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Non-ferroelectric nature of the conductance hysteresis in CHsNHsPbIs 
perovskite-based photovoltaic devices 

J. Beilsten-Edmands/ G. E. Eperon,^ R. D. Johnson/ H. J. Snaith/ and P. G. Radaelli^ 

Department of Physics, Clarendon Laboratory, University of Oxford, Oxford 0X1 3PU, 

United Kingdor^^ 

We present measurements of conductance hysteresis on CHaNHaPbla perovskite thin films, performed using 
the double-wave method, in order to investigate the possibility of a ferroelectric response. A strong frequency 
dependence of the hysteresis is observed in the range of 0.1 Hz to 150 Hz, with a hysteretic charge density in 
excess of lOOOpCcm”^ at frequencies below 0.4 Hz — a behaviour uncharacteristic of a ferroelectric response. 
We show that the observed hysteretic conductance, as well as the presence of a double arc in the impedance 
spectroscopy, can be fully explained by the migration of mobile ions under bias on a timescale of seconds. 
Our measurements place an upper limit of « 1 pC cm“^ on any intrinsic frequency-independent polarisation, 
ruling out ferroelectricity as the main cause of current-voltage hysteresis and providing further evidence of 
the importance of ionic migration in modifying the efficiency of CHsNHsPbla devices. 


Solar cells based on a layer of halide perovskite as 
the active material have recently seen rapid increases 
in efficiency, with certified values in excess of 209^. In 
these devices, a hysteresis effect has been observed in 
the current-voltage (JV) curves, whereby the photocur¬ 
rent depends on the direction of the voltage sweep be¬ 
tween forward bias and short circuit configurations. Un¬ 
derstanding this hysteresis effect is important for char¬ 
acterising the steady-state efficiency of such devices re¬ 
liably under working condition^. Hysteresis has been 
reported to depend on scan rate, contact materiaP, bias¬ 
ing historjBl, and lighting condition^. In addition, time 
dependent photoconductivity responses have also been 
observed^. Several explanations for these effects have 
been proposecP, including charge traps at the perovskite 
surface, ion migratiorP leading to band-bending at the 
interfacd^, and t he int rinsic ferroelectric nature of the 
perovskite materiaP*^*21III]_ 

The possibility of ferroelectricity has been suggested 
since, in general, the perovskite crystal structure is 
known to undergo a range of distortions, including po¬ 
lar distortions, thereby allowing a ferroelectric polari¬ 
sation to develop. This theory has gained momentum 
after a recent crystallographic study of the room tem¬ 
perature tetragonal structure reported the polar Mem 
space groujP^, rather than the previously accepted non¬ 
polar /4/mc?7p^®^. In addition to the potential po- 
larisability of the Pb-I lattice, orientational ordering of 
the methylammonium (MA) molecular dipoles could lead 
to molecular ferroelectricity or antiferroelectricitjB^. Ab 
initio calculations have predicted a ferroelectric polari¬ 
sation (P) of 38pCcm“^ for CHaNHsPbla®, while mea¬ 
surements indicating a f erroelectric response have been 
reported in the literatur^HE^EHiZl However, there has 
yet to be any conclusive evidence of the ferroelectric 
properties of this material. Although not intrinsically 
hysteretic, dielectric effects must also be considered, as 
an applied voltage will induce a dielectric polarisation, 
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which may affect device operation. A dielectric study on 
CHsNHaPbIs revealed a picosecond relaxation process, 
corresponding to dynamic disorder of the MA ionaiSl^ 
whereas other studies have identified a low fr equen cy re¬ 
sponse, also attributed to dielectric relaxatioi J^^ l ^^ l 

We set out to investigate whether ferroelectricity is the 
cause of the observed JV hysteresis in CHaNHaPbla thin 
films. Using the double-wave method (DWMj^, also 
known as the PUND methoepJ^, we observed polarisation- 
electric field (PE) hy steresis loops, in agreement with pre¬ 
vious studieJjJElIIi] However, the giant low-frequency 
magnitude of the hysteretic charge density and its large 
frequency dependence are clear demonstrations that this 
response cannot be due to ferroelectricity. Based on the 
observed charging processes, we attribute the PE hys¬ 
teresis to ionic migration under applied voltage and the 
band-bending that this induces, and we show how this 
process manifests as JV hysteresis. Therefore, we believe 
that most attempted measurements of ferroelectricity re¬ 
ported so far have wrongly labelled the effect observed, 
and that this term should not be applied to the switch- 
able behaviour in methylammonium halide perovskites. 

The samples for hysteresis measurements were pinhole- 
free 350 nm thick CHaNHaPbla thin films sandwiched be¬ 
tween two conductive electrodes. The CHaNHaPbla, fab¬ 
ricated using a technique previously described as vapour- 
assisted solution processin^^, was contacted on one side 
with fluorine-doped tin oxide and on the other with a 
gold electrode. Eull details of fabrication and growth are 
described in the supplementary informat iorPS. The po¬ 
tential ferroelectric response was investigated using the 
DWM, in order to separate out unambiguously hysteretic 
behaviour from conductive and dielectric effects. The 
methodology of the DWM is fully described in Ref. [501 
here we provide a brief explanation of our measurement. 
A voltage pulse was applied across the film, which was 
in series with a large reference capacitor C^ef as shown 
in Eigure [^a). The capacitor integrates the surface 
charge due to ferroelectric and dielectric polarisation {Q 
= P X A, where A is the contact area) plus the conduc¬ 
tive current of the film during the measurement. Gref was 
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FIG. 1. (a) Experimental circuit and (b) applied waveform 
used for the DWM. A voltage amplitude of 2V was used, 
corresponding to an E-field of 5.7kVmm“^. Measurements 
were made across a frequency range of 150 Hz to 0.1 Hz, corre¬ 
sponding to pulse width times tpuise from 1/300 s to 5 s. Gap 
times At between 5 s and 10 s were necessary to discharge Gref, 
(c) Expected E-field dependence of accumulated charges due 
to ferroelectric polarisation, dielectric polarisation and diode 
conduction. For semiconductive materials Qcond can be or¬ 
ders of magnitude larger than Qperro = Pperro x A, where A 
is the contact area. 


chosen to fulfil the critereon i?sampieCref tpuise such 
that additional transient behaviour was not introduced 
(with values ranging from 2.2 mF at 150 Hz to 150 mF 
at 0.1 Hz). This maintained the capacitor voltage below 
0.2 V and this back-voltage was corrected for in the loop 
analysis. The applied voltage waveform is shown in Fig¬ 
ure Bt) and experimental parameters are given in the 
figure caption. The waveform consisted of a sine wave 
of one period, followed by two half period sine waves in 
the positive direction and two in the negative direction, 
hereafter referred to as the initial pulse and pulses 1 to 
4. Gap times ranging between 5 and 10 seconds were 
necessary to discharge Gref between pulses, to avoid the 
introduction of additional hysteretic behaviour. The ini¬ 
tial pulse would pole a ferroelectric film in the ‘down’ 
state. The integrated charge from pulse 1 would contain 
any hysteretic switching behaviour, as well as additional 
contributions from non-hysteretic effects such as conduc¬ 
tion and dielectric polarisation. Figure [^c) shows the 
expected electric-field dependence of these contributions. 
As the him is now in the ‘up’ state, pulse 2 gives only 
the unwanted additional contributions, and subtraction 
yields the underlying hysteretic behaviour. This process 
is repeated for the negative direction, enabling the full 
charge-voltage (QV) hysteresis loop (which can be ex¬ 
pressed as a PE hysteresis loop) to be analytically re¬ 
constructed. All measurements were made under dark 
conditions, to eliminate the photocurrent and measure 
the intrinsic characteristics of the hlms. 

Figure [^a) shows a hysteresis loop obtained by the 
DWM on a CHaNHaPbla thin him. The hysteresis is 


plotted in units of pCcm“^, in order to compare to 
known ferroelectrics. Figures[^b),(c) show the frequency 
dependence of the hysteretic charge density measured un¬ 
der positive bias. Remarkably, the shape of the hysteresis 
loop looks very similar to the one expected for a fer¬ 
roelectric material but the magnitude of the hysteretic 
charge density exceeds lOOOpGcm”^ below 0.4 Hz — 
over an order of magnitude larger than for perovskite fer¬ 
roelectrics. For example, the polarisation of bulk BaTiOa 
is 15pCcm“^ (Ref. [M]), and is enhanced to 70pCcm“^ 
in thin hlm^^, while the largest reported polarisation 
is on the order of 150pCcm“^ for BiFeOa thin hlm^^. 
Hence, our measured hysteresis in GHsNHaPbla is very 
unlikely to arise from ferroelectricity. As explained pre¬ 
viously, this large charge density cannot simply be at¬ 
tributed to ferroelectricity plus an additional contribu¬ 
tion due to electronic conduction, as this is corrected for 
by the DWM. Additionally, between 150 Hz and 4 Hz, the 
measured hysteresis is to a good approximation linear in 
1// (i.e. linear in the time of applied voltage). This is 
uncharacteristic of the switching kinetics observed in fer¬ 
roelectric thin film^^, where ferroelectric polarisation is 
expected to be frequency independent to good approx¬ 
imation, particularly in the low frequency regime. The 
giant magnitude of the hysteretic charge density requires 
a dominant non-ferroelectric process, which, as explained 
in the remainder of the paper, we identify as due to ion 
migration and trapping at the interfaces. Since these 
effects are minimised at high frequency, we can define 
an upper limit to the ‘true’ ferroelectric polarisation as 
the high-frequency intercept of the Qh vs 1// line. As 
1// —>■ 0, the linear fit was found to pass through the 
origin within error, and the uncertainty puts an upper 
limit on any intrinsic frequency-independent polarisation 
of lpCcm“^. Although this uncertainty is not insignif¬ 
icant, we can state with confidence that our measure¬ 
ments rule out a large polarisation as predicted by Frost 
et aZ.^. It should be noted that if there were a ferroelec¬ 
tric molecular dipole ordering that had a slow response to 
an external field (i.e. at low frequencies), this would not 
be visible in our measurements due to the dominant non- 
ferroelectric contribution to hysteresis at low frequency. 
The absence of ferroelectricity is in agreement with the 
majority of reports in the literature on the crystallogra¬ 
phy of GH 3 NH 3 Pbl 3 , which determine a non-ferroelectric 
structure described by the 74/mem space groujP^®^. A 
polar distortion of the Pb-I lattice or ferroelectric order¬ 
ing of MA dipoles would necessarily lower the symmetry 
of the whole crystal to a polar space group. 

In order to establish the true origin of the QV hystere¬ 
sis, the majority of the frequency dependence can be re¬ 
moved by considering the current density, rather than ac¬ 
cumulated charge. Figure [^d) shows the dark JV curve 
as measured by the DWM at 0.2 Hz. The measured hys¬ 
teresis is now visible as a difference between the current 
density extracted from the film on the first and second 
(or third and fourth) pulses, with this JV hysteresis be¬ 
ing frequency independent above 1 Hz, and saturating to 
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FIG. 2. (a) Hysteretic charge density 

(Qh/A) as measured by the DWM on a 
CHgNHaPbls thin film at 2 Hz. (b) Fre¬ 
quency dependence of the QV hysteresis as 
a function of 1// below 4Hz (circles). The 
continuous line is a linear fit to the data for 
1/./ < 0.25. (c) Frequency dependence of 
the QV hysteresis at low frequencies, where 
some saturation from the low frequency lin¬ 
ear behaviour (dotted line) is observed. For 
/ below 1 Hz, error bars indicate the range 
of measured values of hysteresis at each fre¬ 
quency. (d) DWM measurement plotted as 
JV curves. Arrows indicate voltage sweep 
direction. The QV hysteresis arises from a 
difference in current density between pulses 
1 & 2, and pulses 3 & 4. The voltage axis 
has been reversed to display the JV curve 
in the orientation typically shown in solar 
cell literature. 


some extent below 1 Hz. It is this frequency indepen¬ 
dent JV hysteresis above IHz that gives rise to a QV 
hysteresis linear in 1// when the current is integrated in 
the DWM measurement. A reduction in absolute current 
density as the frequency was reduced was also observed. 

To further investigate the origin of JV hysteresis, we 
measured the time-dependent response from the film at 
fixed applied voltage V, after initial poling of -V, repre¬ 
sentative of DC switching. Figure [^a) shows the time- 
dependent current density at a series of applied voltages, 
measured on a single film. We note that the results 
measured on this film, discussed in the following, were 
fully representative of a number of films measured. It 
was not possible to apply constant voltages above 1.5 V 
as this was found to damage the films. In the con¬ 
stant voltage measurements we observed a lower overall 
current density compared to the DWM measurements, 
which we attribute to film degradation over time. At 
fixed voltage, the initially observed current is greater 
than the steady-state current, which is reached after ap¬ 
proximately 10 seconds. This is indicative of a charge¬ 
trapping process at the perovskite/electrode interface 
that saturates on the timescale of seconds. Although this 
could in principle be due to either electronic trapping or 
ionic trap pin g, given the high electronic mobility of this 
compouncp21^ one would expect electronic trapping on 
the order of nanoseconds. Therefore this charge should 
be attributed to ionic migration, which has indeed been 
previously reported to occur on the timescale of seconds 
in CHsNHaPbla devic ej^. P erovskite halides are known 
to be ionic conductor^^HH and the applied maximum 
field of 5.7kVmm“^ should induce significant migration. 
Upon removing the voltage, the discharging current was 
negligible, showing that the additional charge had indeed 
been trapped. To compare directly to our JV and DWM 
measurements, we integrated the additional charge above 


the steady state background at 10 s at each fixed voltage, 
between 0 and 10 s. Figure [^b) shows a comparison of 
the lowest frequency DWM hysteresis loop and a hystere¬ 
sis loop constructed from the integrated excess charge 
after switching at constant voltage. The reduced magni¬ 
tude of accumulated charge in the constant voltage mea¬ 
surement is attributed to current-induced fatigue, as pre¬ 
viously stated. The agreement between the two methods 
reveals how ferroelectric-like hysteresis loops can result 
from slowly saturating ionic charges and ionic trapping at 
the interfaces, where the role of the coercive electric field 
is played by the activation energy for ionic de-trapping. 
Due to the time dependent nature of the voltage wave¬ 
form in the DWM, at later times the ionic current is 
always reduced for a given voltage, resulting in a differ¬ 
ence once the pulses are subtracted. We remark that, 
in addition to charge trapping, any change in the carrier 
extraction efficiency at the perovskite/electrode interface 
due to migration-induced band-bending will modify the 
steady-state background and hence will also give a con¬ 
tribution to the additional charge observed, in both sets 
of measurements. Our measurements also enable us to 
place an upper limit on the fraction of ions required to 
give such an accumulated charge. As an example, tak¬ 
ing the I~ ions to be the ionic species, the hysteresis of 
2500pCcm“^ measured using the DWM corresponds to 
an upper limit of 3.7% of the I~ ions present in the film. 
In reality it is likely that all species give some contribu¬ 
tion to ionic migration since, for CH 3 NH 3 Pbl 3 , vacancies 
of MA, Pb and I have all been calculated to have low for¬ 
mation energies (below 2eV), depending on the chemical 
potentiaP^. It is important to note that whilst our DWM 
measurements start at 0 V, typically JV measurements in 
the literature are started from open circuit, equivalent to 
measuring the second half of pulse 3 followed by pulse 4. 
Therefore this measurement shows an additional feature 
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FIG. 3. (a) Time-dependent cur¬ 

rent density (circles) at a series of 
fixed voltages V, after poling at -V. 
The solid lines indicate the saturated 
steady-state backgrounds, (b) Com¬ 
parison of the integrated additional 
charge at fixed V and DWM mea¬ 
surement at 0.1 Hz. The integrated 
charge is shifted by —2500pCcm“^ 
such that both curves have a com¬ 
mon origin at (0,-2500). (c) Equiv¬ 
alent circuit for our thin films, in¬ 
cluding the effect of ionic migra¬ 
tion. (d) Calculated impedance re¬ 
sponse of the equivalent circuit, us¬ 
ing Ctot = 14 nF and the experimen¬ 
tally determined parameters i?rec = 
24012, Rion = 84012 and Cion = 
0.28 mF. The existence of the ionic 
component in the equivalent circuit 
produces an additional arc at low 
frequencies. 


of JV hysteresis at fast sweep rates due to ionic migra¬ 
tion. We expect to be particularly sensitive to any ionic 
charge density at the perovskite interface due to the ab¬ 
sence of electron/hole transporting layers in our films. 
The observation of JV hysteresis at slow sweep rates in 
the literature could then be due to the fact that over 
long timescales the accumulated ions will migrate back 
through the film, causing a loss of the ionically charged 
interface and associated band-bending. 

In the last part of the paper, we discuss the non- 
hysteretic part of the frequency-dependent conductiv¬ 
ity in CHaNHaPbla, focussing particularly on the low 
frequency behaviour, since previous impedance spec¬ 
troscopy data have shown the presence of a low frequency 
arc centred at 0.28 occurring on a timescale of 

secondJi^. Figure [^c) shows a simple equivalent cir¬ 
cuit for the film, expanding on a previous model^ to 
include the effects of ionic migration whilst neglecting 
any stray capacitance or substrate series resistance. i?rec 
represents the electronic recombination resistance of the 
cell and Ctot is the sum of the chemical and geometric 
capacitances, excluding any ionic effects. The additional 
effect of ionic migration can be modelled by a series RC 
circuit in parallel to the electronic processes. This phys¬ 
ically represents an additional current channel that is 
blocked at the contacts, which, at constant voltage, gives 
an initial current of V/i?ion that reduces exponentially 
to zero as Cion becomes fully saturated. We emphasise 
that simple (linear) RC circuits cannot model the hys¬ 
teresis resulting from tra ppi ng behaviour, which requires 
nonlinear RC componentP^l. From the stored charge den¬ 
sity of 2.5mCcm“^ at 2V (figure [^, one can calculate 
an ionic capacitance Con = 0.28 mF , which is remark¬ 
ably similar to the low frequency capacitance reported 


in Ref. [TH For the parameters i?rec and Rion, we used 
an ohmic approximation to the resistance based on the 
current observed in the DWM at -IV (see figure [^d)), 
giving Rrec = 24012, Rion = 84012. We approximated 
Ctot = 14 nF based on the predicted value of the relative 
permittivity ei.=24 in the literatur^^ and neglecting the 
chemical capacitance. Figure |^d) shows the impedance 
response of the equivalent circuit, calculated without any 
adjustable parameters. We find that the addition of the 
ionic component to the equivalent circuit introduces a low 
frequency arc in the impedance plot, as observed in the 
experimental^, but with a somewhat higher characteristic 
frequency /if = 0.68 Hz. It is important to stress that the 
exact values of the characteristic frequencies are sensitive 
to the geometrical parameters of the film. For example, 
the relative sizes of the intermediate and low frequency 
arcs are determined by the ratio of Rion to Rrec, whilst an 
increase in the magnitude of the resistances will lower the 
characteristic frequencies and increase the width of the 
impedance plot (see figure [^d)). Therefore, we conclude 
that our equivalent circuit model, including the effects of 
ionic migration, provides a good explanation of th e low 
frequency features seen in impedance spectroscopjU^^l^l 
Our model estimates ionic mobilities in this compound 
on the order of 10“^°cm^ V“^ s“^, ionic conductivities on 
the order of 10~^f2“^cm“^ and carrier densities on the 
order of 10^°cm“^. Further studies into the ionic mo¬ 
bilities of different ionic species in this compound would 
be beneficial in confirming this model. We remark that, 
according to our explanation, the low frequency arc is 
purely due to conducta nce ra ther than dielectric effect^^, 
as previously propose d^^b9 | Whilst some dielectric po¬ 
larisation will be induced during the pulse, given the fast 
timescales of MA ion relaxatioiJi^ this can only modify 
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the JV characteristics as a function of voltage and will 
not contribute to the hysteresis observed. 

In summary, we have completely characterised the hys- 
teretic behaviour in CHaNHaPbIs thin films and have 
concluded that it cannot be attributed to ferroelectricity. 
The frequency dependent nature of the hysteretic charge 
is uncharacteristic of ferroelectricity, whilst its magni¬ 
tude at low frequencies is over an order of magnitude 
larger than for well-known perovskite ferroelectrics. Im¬ 
portantly, we have ruled out any ferroelectric response 
larger than 1 pC cm“^, significantly lower than previously 
suggested for this compouncP^. The observed hysteretic 
conduction effects are shown to be the result of ionic mi¬ 
gration and of the associated charge trapping, occurring 
on a timescale of seconds. Consequently, it is not correct 
to label the observed effects in CHaNHaPbla solar cells 
as a ferroelectric response. Understanding the ionically 
induced modifications of the electrode/perovskite inter¬ 
face in this material is likely to be the key to explaining 
the JV characteristics and to enhance the performance 
of perovskite solar cells. 
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Supplementary Information: Non-ferroelectric nature of the conductance 
hysteresis in CHsNHsPbIs perovskite-based photovoltaic devices 
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Materials. Unless otherwise stated, all materials were purchased from Sigma-Aldrich or Alfa Aesar and used as 
received. Spiro-OMeTAD was purchased from Borun Chemicals. Methylammonium Iodide (MAI) was synthesised 
according to a reported procedure.^ 

Substrate Preparation. Devices were fabricated on fluorine-doped tin oxide (FTO) coated glass (Pilkington, 7 D/sq). 
Initially FTO was removed from regions under the anode contact by etching the FTO with 2 M HCl and zinc powder. 
Substrates were then cleaned sequentially in hallmanex detergent, acetone, propan-2-ol and oxygen plasma. 

CH^NH^Pbl^ film formation. Firstly a 0.75 M solution of lead (II) iodide was spin-coated on the FTO substrate at 
3000 rpm for 15 s, with both solution and substrate at 85 °C immediately prior to coating. Films were subsequently 
dried at 85 °C for 10 minutes. Films were then transferred into a nitrogen-filled glovebox and placed on raised supports, 
with the lead iodide facing upwards, in a small sealed glass chamber. 900 mg of MAI was crushed with a pestle and 
mortar and placed on the bottom surface of the chamber. The chamber was then heated from below on a hotplate at 
1555 °C for 5 hours. Films were finally rinsed in anhydrous propan-2-ol and annealed for 30 minutes at 100 °C. 

^M. M. Lee, J. Teuscher, T. Miyasaka, T. N. Murakami, and H. J. Snaith, Science 338, 643 (2012). 



